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We describe a realistic model for a focused high-intensity laser pulse in three dimensions. Rela- 
tivistic dynamics of an electron submitted to such pulse is described by equations of motion with 
ponderomotive potential depending on a single free parameter in the problem, which we refer to as 
the "asymmetry parameter". It is shown that the asymmetry parameter can be chosen to provide 
quantitative agreement of the developed theory with experimental results of Malka et al. [Phys. 
. . . Rev. Lett. 78, 3314 (1997)] who detected angular asymmetry in the spatial pattern of electrons 

■ accelerated in vacuum by a high-intensity laser pulse. 
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^ ' In their recent paper, Malka et al. 0] have reported experimental observation of electrons accelerated to relativistic 
O . energies by a high-intensity linearly polarized subpicosecond laser pulse in vacuum (see also comments on the paper 
' and the author's reply in Ref. (2|-U]). This effect, known as high-intensity ponderomotive scattering, was discussed 

m detail in Ref. §]. It occurs when the quiver amplitude imparted by the laser field to an electron becomes comparable 
to the focal spot radius of the laser beam. If the beam is Gaussian, the radial restoring force acting on the electron 
decays exponentially and the electron can be scattered out of the pulse. 
I ' The data of Malka et al. |^ show that the energies gained by the scattered electrons are in good quantitative 
^ . agreement with calculations of electron trajectories in the polarization plane made with the first-order paraxial model 
' for the laser field Nevertheless, the first order paraxial model predicts isotropic electron scattering [|| that 

is not supported by experimental results. Indeed, accelerated electrons were detected by Malka et al. only in the 
(E,k) plane, while no significant signal was detected after rotating the laser polarization direction by 90° In 
our opinion, this discrepancy between the theory and experiment is due to the following. In the first-order paraxial 
model focusing of a plane monochromatic wave only leads to an appearance of nonvanishing longitudinal components 
of electromagnetic fields in the focal region. However, the focusing is known to affect transverse components of the 
^ ' fields also (see, e.g. Ref. [^). In particular a plane monochromatic wave, polarized linearly along the x axis and 
Oh, propagating along the z axis, is converted by an aplantic system to a converging spherical wave with non-vanishing y 

and X components of electric and magnetic fields, respectively. 
^ , In this Letter, we show that experimentally observed 0] anisotropy of ponderomotive electron scattering can be 
explained in the framework of a realistic model for the laser field developed in our recent paper Q. Our model is 
based on an exact solution of Maxwell equations in three dimensions (3D), which can serve to describe a stationary, 
focused monochromatic laser beam with characteristic frequency to and arbitrary intensity. Amplitudes of electric and 
magnetic fields in the model depend on radial coordinates as well as the coordinate along the direction of the beam 
propagation. These amplitudes are characterized by parameters R and L = ujR^, which can be interpreted as the focal 
spot radius and the Rayleigh length of the laser beam, respectively. The model admits different field configurations, 
which are determined by two coordinate functions satisfying certain second-order partial differential equations. Some 
special choice of these functions describes the Gaussian beams, which are widely used in optics. The model can 
be generalized by introducing temporal amplitude envelope g{ip/u)T), where if is the relativistically invariant phase 
of the traveling wave, to describe a laser pulse with finite duration r. (It is assumed that the function g(ip/ujT) is 
equal to unity at the point ip — and decreases exponentially at the periphery of the pulse for \ip\ ^ ujt.) In this 
case the electric and magnetic fields of the model constitute an approximate solution of Maxwell equations with the 
second-order accuracy with respect to small parameters A and A', defined as 

A' = 1/wT < A = l/tji?< 1. (1) 
For a pulse propagating along the z direction, x and y components of the electric field oscillate with phase difference, 
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which depends on z and values of x and y coordinates of a point in the plain z = const. Moreover, the aforementioned 
phase difference depends also on ip. Therefore, one cannot ascribe some definite type of polarization to a tightly 
focused laser pulse. Nevertheless, for a weakly focused pulse (A <C 1), there always exists a region near the axis of the 
beam r R, where the field properties are very close to those of a plane wave field. This region we call "the plane 
wave zone" . It is reasonable to ascribe polarization of the field in this region to the beam as a whole. Hereafter we 
refer to the field of the pulse as linearly polarized in this sense only. For a tightly focused beam the focal spot radius 
is of the order of the wavelength and the plane wave zone doesn't exist. Therefore, only polarization of the parental 
beam incident on the focusing optical system can be ascribed to the focused beam in this case. 

An arbitrary field, linearly polarized along the x axis, may be represented |^] as a superposition of E- and H- 
polarized waves (i.e. waves with the vectors E and H being perpendicular to the direction of the pulse propagation, 
compare with Ref. ||] ) . Relative contributions of E- and i?-polarized waves to the resulting field are characterized by 
the "asymmetry parameter" jj, 



, S" ' -oo<M<oo, (2) 

where E^'q are the x components of the electric field for E- and -ff-polarized waves at the focal point r = for iy9 = 0. 
Note, that in contrast to the amplitude, the quantities E^'^ can take both positive and negative values. 

In the lowest approximation in A and A', the averaged equations of motion of electrons (ponderomotive equations) 
in the field of the linearly polarized laser pulse take the form j^] 
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Here — (p^), where p'' is the 4- momentum of the electron, q^ = qo — qz, p± = {r/R), ( ~ (z/R), brackets () mean 
averaging over fast oscillations, and the ponderomotive potential U is defined by the expression 

U = ^^gHv/^r) { |Fi I V |i^2 1' + M cos 2^ {F,F* + F*E,) } , (4) 
where tanip = Qxl Qy, and ryo is the value of the dimensionless field intensity parameter 

at the focal point at the moment ip — Q {r\Q — a/\/2, where a is the parameter of Malka et at). Functions 
Fi{g±, C; A) are chosen in the form corresponding to the Gaussian beam Q 
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The equation shows that the ponderomotive potential U depends on the azimuthal angle ip, and hence is 
generally speaking asymmetric. The potential U is symmetric only for the case /i = 0. The shape of the ponderomotive 
potential in the plane z = for = is shown in Fig.^ for the cases /i = and p = —1.55. Figure represents 
the ponderomotive potential for the standard case of Gaussian beam commonly used in literature, while Figure ^& 
illustrates the dramatic difference between the cases p = and p ^ 0. For p Q, the ponderomotive potential, 
possesses (besides the central peak) two extra maxima, which are located in the polarization plane. They arise as a 
result of the non-uniform intensity distribution in the plane z = for p ^ Q. Locations of the additional maxima, as 
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well as their amplitudes, are determined by the value of ^. It is noteworthy, that the case ^ = is the only one when 
Ey- and iJa;-components of the electric and magnetic fields remain to be equal to zero outside the plane wave zone 
for the pulse polarized along the x axis. 

We use Eqs. and (Q) for our analysis of free electron acceleration by a co-propagating intense laser pulse in 
vacuum under conditions close to those used in the experiments of Malka et al. [Q . The initial electron energy is taken 
to be e = 10 keV (wo — 0.2c). The laser field parameters are: A = 1 /im, 770 — 2.12 (corresponds to the parameter 
a = 3 of Malka et ai), R = 10 /j,m and ojt = 480. For the temporal envelope of the pulse g{(p/u!T) a sine-squared 
shape is taken. The asymmetry parameter is not determined experimentally and remains a free parameter of the 
problem. Its value, fi — —1.55, has been chosen for better fitting of our computational results to the results of the 
experiment. 

It is clear that the maximum energy will be gained by electrons that initially propagate exactly along the axis of 
the laser beam. However, the ponderomotive equations yield zero net energy transfer for such electrons, since they 
can "feel" the spatial gradient of the Gaussian laser field, and hence can be scattered out of the pulse, only due to 
the quiver motion which is absent for the average trajectory described by the ponderomotive equations. Nevertheless, 
there exist a family of trajectories with nonzero initial distances from the beam axis, for which the gained energy is 
close to its maximum value. The energy and scattering angle of the electrons depend also on the position at which 
the particle is overtaken by the pulse [Q . Indeed, the maximum energy is obtained for the electrons that experience 
the peak field of the laser and therefore meet the laser pulse at some distance before the focus. As a result, we obtain 
a 3D domain of injection positions of electrons obtaining the final kinetic energy W > 0.9 MeV at the scattering angle 
39.5°. 

Different cross-sections of this domain for the potential with fi = —1.55 (Fig. |l]&) are shown in Fig. ||. The 
longitudinal size of the domain is of the order of the Rayleigh length for the laser beam L, whereas its transverse 
size is much less than the focal spot radius R. The cross sections in Fig. || display high degree of radial anisotropy. 
Their shape, of course, is essentially determined by the type of the ponderomotive potential or by the value of jj,. 
In particular, the domain of injection positions for the case (not presented here for the sake of compactness) of the 
potential shown in Fig. with /x = is purely radial, in agreement with 

To obtain the angular distribution of scattered electrons we have calculated their trajectories numerically, applying 
the fourth-order Runge-Kutta method to the ponderomotive equations (|^). Initial positions of the electrons were 
taken from the domain shown in Fig. ^. Since the shape of the cross-sections z — const varies very slowly at intervals 
Sz ^ R, we have considered equidistant planes z = nR with n = —27, —26, . . . , 5. In each of these cross-sections, 
the electron injection positions were chosen randomly under condition that their density was constant and equal to 
3 X 10^^/i?^. Physically, such procedure corresponds to uniformity of the initial electron beam. The total number of 
electron trajectories considered in such a way was more than 2 x 10^. We were interested only in those trajectories, 
which crossed the plane z = 11.66 cm at the points contained inside the ring with radii ri = 8.99 cm and r2 = 9.89 cm. 
The latter conditions were determined by the position and angular size of the detector in the experiment 0] . 

The results of the calculations are presented in Fig. ^. We plot the normalized number of scattered electrons (n) 
with final energies W > 0.9 MeV as a function of the azimuthal scattering angle a. The values of (n) at any given a 
were obtained by averaging the number of scattered electrons over the range of a equal to the angular size of the 
detector used in the experiments 0] . It is easily seen, that the angular distribution of scattered electrons essentially 
depends on the parameter of asymmetry /i. At fi — 0, which corresponds to the symmetrical Gaussian ponderomotive 
potential shown in Fig. |l|a, the distribution is isotropic and purely radial (compare to the result of Ref . ^ ) . At the 
same time, if /i = —1.55, the number of scattered electrons detected in the (E,k) plane {a = 0) is about 30 times 
higher than that in the (H,k) plane (a — 7r/2). This result is in good quantitative agreement with observations of 
Malka et al. 0. 

The 30-fold anisotropy of accelerated electrons is clearly explained by asymmetry of the ponderomotive potential 
. The cross term in (Q), besides asymmetric corrections to the radial force, gives rise to a tangential force which 
is responsible for pushing electrons out of the plain perpendicular to the polarization plane. The latter corresponds 
to minimum of the ponderomotive potential as function of azimuthal angle ip, while the perpendicular plane to its 
maximum. Therefore one could be surprised that the number of electrons scattered at the angle a = 7r/2 is not equal 
to zero. Certainly it is explained by complicated structure of the ponderomotive potential (^, namely by the fact 
that the cross term can change its sign at the periphery of the focus. 

The asymmetry of the ponderomotive potential itself is determined by non-zero value of the parameter /i which 
characterizes relative contributions of E- and H-polarized waves to the resulting field. As far as we know, nobody 
has never controlled the parameter fi in experiments. The reason is evident. Before the work of Malka et al Q there 
were no experiments where the three-dimensional intensity distribution of a laser pulse influenced physical results. 
Therefore even qualitative coincidence of our calculations with the results of the experiment M would give a deeper 
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insight into physics of the electron-laser interaction. However, it appeared that the model describes the experiment 
quantitatively. Of course, the quantitative agreement between the developed theory and the experiment is based on 
fitting of a single free parameter fj, in the problem, which has not been measured experimentally. Therefore, from the 
standpoint of our model, the experiment of Malka et al. |^ could be considered as a probe for 3D field distribution in 
the laser pulse. The correctness of our approach could be verified by an independent experiment for another physical 
situation performed with the same laser system. Measurements of angular distribution of ATI electrons could serve 
as a good example of such experiment. 

We thank M.V. Fedorov and V.D. Mur for fruitful discussions. This work was supported by the Russian Foundation 
for Basic Research under projects 00-02-16354 and 00-02-17078. 
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FIG. 1. Ponderomotive potential (^ in the plane z = at ^ = for two values of asymmetry parameter:(a) ^ — Q and (b) 
= -1.55. 

FIG. 2. Cross sections of 3D injection-position domain for electrons gaining the energy W > 0.9 MeV in the ponderomotive 
potential (^ with /j, — —1.55. z coordinates of the cross sections are (a) z = — 21ii, (b) z — —167?, (c) z = — 11_R, (d) z = — 6-R, 
(e) z = ~R, and (f) z = 2R. 



FIG. 3. Normalized number of scattered electrons with the energy W > 0.9 MeV as a function of azimuthal angle a for (a) 
fi — and (b) /i — —1.55. 
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